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By using topological current theory, we study the inner topological structure of disclinations 
during the melting of two-dimensional systems. From two-dimensional elasticity theory, it is found 
topological currents for topological defects in homogeneous equation. The evolution of disclinations 
is studied, and the branch conditions for generating, annihilating, crossing, splitting and merging of 
disclinations are given. 
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Topological defects, which are a necessary consequence 
of broken continuous symmetry, play an important role 
in two-dimensional phase transition. In 1970's, Koster- 
litz and Thouless construct a detailed and complete the- 
ory of superfluidity on two-dimensions [l|. They in- 
dicate vortices pair unbinding will lead to a second- 
order transition in superfiuid films. Later, a micro- 
scopic scenario of 2D melting has been posited in the 
form of the Kostcrlitz-Thouless-Halperin-Nelson- Young 
(KTHNY) theory @,i,H. The KTHNY theory predicts 
a new phase, the so-called hexatic phase, that exists be- 
tween the solid and liquid phases in 2D melting 

In two-dimensional colloid systems, topological defects 
have been studied in experiments and computer simu- 
lations. A serial experiments were performed to cal- 
culate dislocations and disclinations dynamic of two- 
dimensional colloidal systems, and dissociation of dislo- 
cations and disclinations were observed 

0,0,11, EH. D^- 
ing the years, a large number of computer simulations 

indicated that exist a two-stage melting as prescribed 
by KTHNY theory, however, results are still controver- 
sial [HE1E1E1- Althou S h the KTHNY theory is cur- 
rently preferred, a different theoretical approach, evoking 
grain-boundary-induced melting, was a first-order tran- 
sition suggested by Chui [lj|. One may note that the 
condensation of g eometrical defects is also a first-order 
transition [l7l Il8j |. Our previous work found that exist 
a hexatic-isotropic liquid phase coexistence during the 
melting of soft Yukawa systems fl9l . [20| . By Voronoi 
polygons analysis, the behavior of piont defects in the 
coexistence is very complicated. The evolution of topo- 
logical defects during the melting of two-dimensional sys- 
tem still a open question. 

Recently, a topological field theory for topological de- 
fects developed by Duan et al 21]. By using Duan's topo- 
logical current theory, the inner topological structure and 
bifurcation of topological defects, such as disclination and 
dislocation in liquid crystal and soild, were studied. In 



KT phase transition, there also exists an elementary vor- 
tex topological current constructed by the superfiuid or- 
der parameter [22|]. By using the topological current the- 
ory, we can give the the branch conditions for generating, 
annihilating, crossing, splitting and merging of topologi- 
cal defects. 

In this paper, we will discuss the topological quan- 
tization and bifurcation of topological defects in two- 
dimensional crystals. This work is based on the so-called 
Duan's topological current theory. The organization of 
this paper is as follows. We describe the elasticity theory. 
Using Duan's topological current theory, we discussed the 
topological structure of disclination in two-dimensional 
crystals. In the last section, we summarize our results. 

In continuum elasticity theory, the elastic Hamiltonian 
in two-dimensional triangular solid is given by [23| 



F 



d 2 r(2fiu 2 lj + Xu 2 kk ), 



(1) 



where A and /j. are the two-dimensional Lame coefficients. 
The strain tensor is 



u ij( r ) = o 



dui(r) duj(r) 
drj dri 



(2) 



A deformation is represented by a displacement vector 
field u(r) = {u\,u<z), which maps the point r = (x,y) 
to r + u. If there are no defects, the deformation is a 
single-valued mapping of the plane onto itself. But u 
becomes a multi-valued function when there is a dislo- 
cation. A single dislocation corresponds to an extra half 
lattice plane, which characterized by a Burger's vector b. 
Another type of defect in two-dimension crystal is discli- 
nation, which is defined in terms of the bond angle field 
6. 6{x) is the angle between local lattice bonds and a 
reference axis. 

If we minimize F s with respect to variations in u, we 
obtain the equation 
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di(Jij = 0, 

Where the stress tensor <jij is defined by 
= 2fiUij + \u kk 5ij, 



(3) 
(4) 
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Because er^ is symmetric, it can be written as 
<Jij = eikCjiVk^iX, 



(5) 



The function \ is called the Airy stress function. Al- 
though any choice for \ yields a stress tensor that sat- 
isfies Eq.([3]), the choice cannot be arbitrary. The strain 



Uij is related to the stress is 



1 



A 



Sin cr kk 



uy " 2/T y 4(A + /i) 11 

_ 1 e ik e kl V k V ix - yV 2 X(5 y (6) 

where Y = 4/i(^i + A)/(2/i + A) is the two-dimensional 
Young's modulus and cr = A/(2^t + A) the two- 
dimensional Poisson ratio. Applying e^ej/VfeV; to both 
of this equation, we find 



e k id k di8 + eikdkiejididjU,, 



(7) 



The defects associated with the continuum elastic the- 
ory of a solid are dislocations and disclinations. Disloca- 
tions and disclinations can be introduced into the theory 
in a way similar to the discussion of superfluid vortices [? 
] . In the following we consider only the triangular lattice 
since it is the most densely packed one in two-dimensional 
and favored by Nature. 

Disclinations, which are characterized by a topologi- 
cal charge, have a much higher energy than dislocations. 
They arc defined in terms of the bond angle field 9(r), 
which measures the bonds orientation. It is convenient to 
define an order parameter for bond orientations, which 
for the triangular lattices is 4>{r) = ipoe t6e ( r \ However, 
the bond angle field is undefined at the disclination cores, 
i.e., the zero points of the order parameter. We rewrite 
the orientation order parameter tp(r) = 4>\ + i<j)\ instead 
of ij){r) = ipo^ 66 ^ ■ Let us define the unit vector field n 
as 

, \\<k\\ = VKK, a = 1,2. (8) 
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Obviously, n a n a = 1. The topological defect is related to 
the zero points of the two-component vector parameter 
*, i.e., 



ct>l(x,y) = 0, 4>l(x,y) = 0. 



(9) 



Suppose there is a defect located at Zi, the topologi- 
cal charge of the defect is defined by the Gauss map n: 



W((j>6,Zi) = ^- <£ e ab n a dn b 
2n ^Ei 



(10) 




FIG. 1: Disclination in triangular lattices. A — n/3 disclina- 
tion with its sevenfold coordinated site in the center (left). A 
7r/3 disclination with its fivefold coordinated site (right). 



Using the stokes' theorem in the exterior differential 
form, one can deduce that 

W(<f> e ,Zi) = ±- I e ab e ij d i n a d j n b d 2 x (11) 



2tt 



We can deduce a topological current of disclinations in 
two-dimensional crystal, 

0Lc = l^taAn^nl = S 2 (cf 6 )J k (^) (12) 

It is the ^-mapping current for disclination. where 
J k {(f)^/x) is the vector Jacobians of 4>e, 



(13) 



It is easy to see that this topological current is identi- 
cally conserved, i.e., 



d k j k = 0. 



(14) 



According to the implicit function theorem, if Jacobian 
determinant 



X X 



(15) 



the solutions of Eq.© can be generally expressed as 

x = xi{t), y = yi(t), l = l,2,...,N, (16) 

which represent N zero points zl(t) (1=1,2,. ..,N) or world 
line of N diclinations D[ in space-time. 

With the (^-function theory, 8 2 {<j)) can be expanded as 



N 



S 2 ( 



o = E 



— ^ \J{<Pq/x) 2 



rS'ir-m) (17) 



where the positive integer Pi is called the Hopf index of 
map x — > <j). The meaning of Pi is that when the point 
r covers the neighborhood of the zero zl once, the vec- 
tor field 06 covers the corresponding region for Pi times. 
Using the implicit function theorem and the definition of 
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FIG. 2: Generating and annihilating of disclination pairs. 



vector Jacobians (|13j) , we can get the velocity of the 1-th 
defect, 



v, 



dzj _ r J((f) 6 /x) -, 
dt ~ hifo/x)** 



J(^) = [J 1 (<fie/x),J 2 (<fie/x)] 



(18) 



Then the spatial and temporal components of the defect 
current j u , can be written as the form of the current and 
the density of the system of N classical points particles 
with topological charge Wi = (3ii]i moving in the (2+1)- 
dimensional space-time, 



N 



N 



i=i 



where r\i is Brouwer degree, 

J((j>e/x) 



m 



\J(fa/x)\ 



±1 



(19) 



(20) 



(21) 



For disclinations, using Duan's topological current the- 
ory the homogeneous equation can write as 



V 4 x 



2tt 
— < 
6 



(22) 



Similar results get by Nelson in the KTHNY theory [24}. 
In our theory, the topological charge of a disclination Di 




FIG. 3: Two disclinations collide with different directions of 
motion at the bifurcation point in (2+l)-dimensional space- 
time. 



is 



Qi = If e ab e^d l n a 6 d ] n\d 2 x 



2tt 



Wi = ftm 



(23) 



where Wi is the winding number of \& around D;, the 
above expression reveals distinctly that the topological 
charge of disclination is not only the winding number, but 
also expressed by the Hopf indices and Brouwer degrees. 
The topological inner structure showed in Eq. (|23|) is more 
essential than that in Eq. ([7]) , this is the advantage of our 
topological description of the disclination. 

It is clearly seen that Eq. (fl9|) shows the movement of 
two-dimension crystal topological defects in space-time. 
According to Eq. (fT4jl , the topological charge of defects in 
two-dimensional crystal are conserved, 



dp 
dt 



Vj = 0. (24) 
In addition, there is a constraint of " charge neutrality" , 



1 N 

pd 2 x = -J2Pim = o 
J i=i 



(25) 



It indicate that the defect in two-dimensional crystal ap- 
pear in pair. 

Analogy to vortices in superfluid films, the zero point 
of the order parameter field play an important role in de- 
scribing the topological defects in two-dimensional crys- 
tal. Now we study the properties of these zero points. If 
the Jacobian determinant J°(4>e/x) ^ 0, we will have the 
isolated solution ofthe zeros of the order parameter field. 
But when J (4>q/x) = 0, the above results will change in 
some way, and will lead to the branch process of defects. 
We denote one of the vectors Jacobian at zero points as 
(t*,zi). According to the values of the vector Jacobian 
at zero points of the order parameter, there are limit 
points and bifurcation points. Each kind corresponds to 
different cases of branch processes. 
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FIG. 4: Two disclinations collide with different directions of 
motion at the bifurcation point in (2+l)-dimensional space- 
time. 



Let us explore what will happen to the disclination at 
the limit point (t*,z{). The limit points are determined 
by 



(26) 



(27) 



Considering the condition (|26|) and making use of the 
implicit function theorem, the solution of Eq.([9]) in the 
neighborhood of the point (t*,zi), 



t = t(x), y = y(x) 
where t* = t(zf). In this case, one can see that 



dx 
~dt 



CO, 



(28) 



(29) 



dt 
dx 



(30) 



The Taylor expansion of t = t(x) at the limit points 

(t*,zl) is 



t-t* 



idH 
2dx^ 



(31) 



which is a parabola in the x-t plane. From this equation, 
we can obtain two solutions x\(f) and X2(t), which give 
two branch solutions (World lines of disclinations). If 



cPt_ 

dx 2 



>0, 



we have the branch solutions for t > t* . It is related 
to the origin of a disclination pair. Otherwise, we have 
the branch solutions for t < t*, which related to the 
annihilation of a disclination pair. 

Since the topological current is identically conserved, 
the topological charges of these two generated or annihi- 
lated disclinations must be opposite at the limit points, 
say 



(hm + fcm = o. 



(32) 



For a limit point it is required that J 1 (0/x)|t* iZ - 7^ 0. 
As to bifurcation point, it must satisfy a more complex 
condition. This case will be discussed in the following. 

Now, let us turn to consider in which the restrictions 
on zero point {t* ,z{) are 



x 



0, k — 0,1, 2, 



(33) 



which imply an important fact that the function relation- 
ship between t and x or y is not unique in the neighbor- 
hood of the bifurcation point (t*,Zi). This fact is easily 
seen from 



dx J 1 {4>q/x) 
dt J(4>q/x) 



dy _ J 2 {(f> 6 /x) 
dt J{4>q/x) 



(34) 



which under Eq. ([33|) directly shows the indefiniteness of 
the direction of integral curve of Eq. (f3~4"]) at (t* , z\). This 
is why the very point (t* , zl) is called a bifurcation point 
of the orientation order parameter. 

With the aim of finding the different directions of all 
branch curves at the bifurcation point, we suppose 



dy 



?0. 



(35) 



(*-vi) 



According to the (/)-mapping theory, the Taylor expansion 
of the solution of the zeros of the order parameter field 
in the neighborhood of (t* ,z{) can be expressed as 

A(x- z\f + 2B(x- z\){t-t*) + C{t-t*) 2 + ... = 0, (36) 
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which includes three important cases, (i) Two world lines 
tangentially contact, i.e., two disclinations tangentially 
encounter at the bifurcation point, (ii) Two world lines 
merge into one world line, i.e., two disclinations merge 
into one disclination at the bifurcation point, (iii) One 
world line resolves into two world lines, i.e., one disclina- 
tions splits into two disclinations at the bifurcation point. 

Case 3(i = 0,C/= 0) For A = A(B 2 - AC) = 
from Eq. ([3"Tj) we have 



dt 
dx 



-B ± VB 2 - AC _ 2B 

~-°> ~~c~- 



C 



(41) 



There are two important cases: (i) One world line re- 
solves into three world lines, i.e., one disclination split 
into three disclinations at the bifurcation point, (ii) 
Three world line merge into one world line, i.e., three 
disclinations merge into one disclination at the bifurca- 
tion point. 

Case 4 (A=C=0). Equations ([37]) and (J3SJ) give re- 
spectively 



FIG. 5: Two disclinations collide with different directions of 
motion at the bifurcation point in (2+l)-dimensional space- 
time. 



which leads to 



and 



(37) 



(38) 



where A, B, and C are constants determined by the or- 
der parameter. The solutions of Eq.([3"T|) or Eq. ([3"8|) give 
different directions of the branch curves (world line of 
vortices) at the bifurcation point. There are four possi- 
ble cases, which will show the physical meanings of the 
bifurcation points. 

Case 1 (4 / 0). For A = A{B 2 - AC) > 0, from 
Eq. (|3T|) we get two different motion directions of the core 
of disclination 



dx I 
~dt I 1 - 2 



-B±^/B 2 ~~AC 

A ' 



(39) 



where two world lines of two disclination intersect which 
different directions at the bifurcation point. This shows 
that two disclinations encounter and the depart at the 
bifurcation point. 

Case 2 (4 / 0). For A = 4(B 2 - AC) > 0, form 
Eq. ([37|) . we obtain only one motion direction of the core 
of disclination 



dx 
~dt 



1,2 



B 

A' 



(40) 



dx 
~db 



o, 



dt 
dx 



(42) 



This case shows that two worldlines intersect normally at 
the bifurcation point, which is similar to case 3. It is no 
surprise that both parts of Eq. (|4"2"[) are correct because 
they give the slope coefficients of two different curves at 
the same point (t*, ij). 

In conclusion, we study the inner topological structure 
of disclinations in two-dimensional colloid systems. We 
have obtained a more essential topological formulary of 
charge density of disclinations in two-dimensional crys- 
tals, and revealed the inner topological relationship of the 
charge of disclinations is characterized by the Hopf index 
and the Brouwer degree. We have studied the evolution 
of disclinations by making use of Duan's topological cur- 
rent theory. We concluded that there exist crucial cases 
of branch processes in the evolution of disclinations when 
J((/)/x) — 0, i.e, rji is indefinite. It means that discli- 
nations are generated or annihilated at the limit points 
and are encountered, split, or merge at the bifurcation 
points, which shows that the disclination is unstable at 
these branch points. We would like to pointed that all 
the results in this paper obtains from the viewpoint of 
topology without any hypothesis, and they are not de- 
pended on the property of systems, such as interaction 
between particles. 
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